Dissipative Hydrodynamic Effects on Baryon Stopping 
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The quark-gluon plasma is considered to behave as a relativistic viscous fluid in the high-energy 
heavy ion collisions. In this study, I develop and estimate a second order dissipative hydrodynamic 
model at finite baryon density with effects of baryon dissipation together with those of shear and 
bulk viscosities. It is found that the hydrodynamic evolution effectively reduces baryon stopping, 
suggesting that the collisions are less transparent at the initial stage. Also the net baryon distribution 
is found sensitive to baryon dissipation as well as to viscosities. The results indicate that the 
dissipative hydrodynamic modeling would be important for understanding unique properties of the 
hot medium even in the high-energy collisions. 

PACS numbers: 25.75.-q, 25.75.Nq, 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

The determination of the properties of the quantum 
chromodynamics (QCD) matter in wider temperature 
and chemical potential regions has been one of the most 
important goals in the hadron physics. The quark mat- 
ter in high temperature with low baryon chemical poten- 
tial is experimentally accessible through the high-energy 
heavy ion collisions, which contributes significantly to 
the exploration of the QCD phase diagram. The early 
experiments predate to the Alternating Gradient Syn- 
chrotron (AGS) in the Brookhaven National Laboratory 
(BNL) and the Super Proton Synchrotron (SPS) at the 
European Organization for Nuclear Research (CERN). A 
sufficiently large amount of energy for the production of 
the quark-gluon plasma (QGP) [l|, a deconfined state of 
quarks and gluons, is considered to be available in the 
^Jsnn = 200 GeV Au-Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) in the BNL , which is one of 
the biggest achievements in the hadron physics. With the 
beginning of the highly-anticipated ^snn = 2.76 TeV 
Pb-Pb collisions at the Large Hadron Collider (LHC) at 
the CERN, the heavy ion programs continue to explore 
the high energy frontiers of the QCD matter. 

At the mid-high-energy collisions in the AGS, the SPS, 
and the early RHIC, the baryon stopping has been used 
to quantify the nuclear transparency in the collisions 
and the kinetic energy loss for the production of the hot 
medium. It is found that the average rapidity losses [!, H| 
at the RHIC are apparently less than the simple linear 
scaling from the AGS (fl and the SPS @ above around 
= 62 A GeV 0,H, which suggests that the col- 
lisions would become increasingly transparent with the 
collision energy. The net baryon distribution itself would 
be an important observable in the higher energy colli- 
sions at the late RHIC and the LHC as shown in a num- 
ber of theoretical analyses even though the net 
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baryon is often treated as irrelevant and neglected be- 
cause it becomes small near mid-rapidity as the collision 
energy increases. It is worth-mentioning that the net 
baryon number in the initial colliding nuclei is conserved 
throughout the time evolution and does not vanish from 
the system even at highest energies. The net baryon dis- 
tribution might also carry valuable information on the 
valence quarks in the colliding nuclei and on the forma- 
tion of the QGP itself in the early thermalization stage, 
which is not well understood yet. 

It i s q uantitatively shown in the RHIC [2(| and the 
LHC [21] experiments that the relativistic hydrodynamic 
analysis is a powerful method to describe the dynami- 
cal behavior of the hot medium when its temperature 
is around and above the crossover temperature [22T - [2~ij ] . 
The typical timescale for the hydrodynamic applicabil- 
ity is t 1-10 fm/c. Recent progresses in hydrody- 
namic studies [25l427j | incorporate non-equilibrium effects 
of shear and/or bulk viscosity into their calculations to 
quantitatively understand the properties of the QGP to 
explain particle spectra, azimuthal anisotropy, multiplic- 
ity, and rapidity distribution. On the other hand, the 
finite density effects are neglected in most of the modern 
hydrodynamic calculations with viscosity although they 
could be important in the context of precision physics. 
It is one of its great advantages that relativistic hydro- 
dynamics can handle finite chemical potentials with ease 
even in off-equilibrium systems. The comparison of the 
experimental data to the hydrodynamical calculations of 
the net baryon distribution would, therefore, provide us 
with valuable insights for understanding the early dy- 
namics as well as for determining the equation of state 
and the transport coefficients, such as baryon diffusion 
coefficient, at non-vanishing baryon density. The devel- 
opment of a consistent dissipative hydrodynamic model 
at finite density would also help the efforts to explore the 
QCD phase diagram in the search for the critical point. 

In this paper, I would like to estimate the effects of col- 
lective flow on the net baryon distribution and the aver- 
age rapidity loss for the late RHIC and the LHC includ- 
ing the interplay of shear viscosity, bulk viscosity, and 
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baryon dissipation. A relativistic dissipative hydrody- 
namic model for the longitudinal expansion of the quark 
matter is developed from a generalized version of the sec- 
ond order theory [28j . Here viscous hydrodynamics refers 
to the hydrodynamics with shear and bulk viscosities, 
which are the tensor and the scalar off-equilibrium pro- 
cesses, and dissipative hydrodynamics to the one with 
charge dissipations, the vector processes, in addition to 
the viscosites [2!|. The initial energy and the initial net 
baryon distributions are employed from the color glass 
theory [30M32l | and the equation of state at finite baryon 
density is constructed from lattice QCD results. 

The paper is organized as follows; Sec. HI is devoted for 
the formulation of relativistic dissipative hydrodynamics 
with finite net baryon density. In Sec. IIII1 the models 
for the equation of state and the transport coefficients 
for hydrodynamic calculations are introduced, along with 
initial conditions and freezeout. In Sec. lIVl I calculate ra- 
pidity distributions of the net baryon number to discuss 
the hydrodynamic effects with baryon dissipation as well 
as shear and bulk viscosities on the baryon stopping. Fi- 
nally in Sec.|Vl discussion and conclusions are presented. 
The natural unit c = ft=/ce = l and the Minkowski 
metric g^ v = diag(+, — , — , — ) is used throughout this 
paper. 



II. RELATIVISTIC DISSIPATIVE 
HYDRODYNAMICS 

I develop a finite-density relativistic hydrodynamic 
model with off-equilibrium processes to describe the dy- 
namical evolution of the quark-gluon plasma in high- 
energy heavy ion collisions. Baryon dissipation is taken 
into account for evolving flow for the first time, together 
with shear viscosity and bulk viscosity. A generalized ver- 
sion [29| of the second order dissipative hydrodynamics 
[28j is introduced for the equations of motion for the dissi- 
pative currents to preserve causality and stability. Cross 
terms among different thermodynamic forces are present 
in the formalism. I focus on the longitudinal evolution 
of the tar get systems and integrate out the transverse 
dynamics [33l l34j in the study because the net baryon 
distribution is expected to depend mainly on rapidity as 
it is carried by the remnant of the valence quarks in the 
colliding nuclei at forward rapidity. It is also experimen- 
tally supported that the net baryon distributions do not 
depend on transverse geometry [35[ . It should be em- 
phasized that the longitudinal boost invariance [361 ] is 
not assumed. 

The energy-momentum tensor T^ v and the net baryon 
current are introduced as the conserving quantities 
in the system. In principle one can introduce other con- 
served currents such as the net strangeness current Njf, 
but those currents are expected to be much smaller than 
Ng for the standard nucleus-nucleus collisions and are 
not considered here. The conserved quantities are re- 
lated to the thermodynamic quantities through tensor 



decomposition as 

T^ v = (e + 8e)u^u v - (P + II) 

+ 2W^U V) + 7T^ , (1) 

N£ = (n B0 + 5n B )u» + V», (2) 

where i// 1 is the flow and A^ = — u^u v is the projec- 
tion operator perpendicular to the flow, i.e., A^u^ = 0. 
In the local rest frame, eo, Pq, and Ubo are interpreted as 
energy density, hydrostatic pressure, and baryon number 
density. II is bulk pressure, energy dissipation cur- 
rent, ir flL ' shear stress tensor, and V* 1 baryon dissipation 
current. The distortions of energy and baryon number 
densities be and bnso equal to zero because of the re- 
quirement of thermodynamic stability (37j . 

The dissipative corrections in the linear response the- 
ory can be interpreted rather intuitively. The bulk pres- 
sure is a dynamical correction to the hydrostatic pressure 
which arises when the system is expanded/compressed 
without deformation, or when the temperature T or the 
chemical potential [Ib decreases/increases. It is note- 
worthy that cancellation among the linear terms could 
be the reason for general smallness of the bulk pressure 
since the phenomena do not occur independently in hy- 
drodynamic systems. See Appendix [A] for the details. 
The shear stress tensor corresponds to the response to 
the deformation without volume change. The energy and 
the baryon dissipation currents are the local fluxes of en- 
ergy and baryon densities which dissipate away from the 
flow, respectively. They are induced by the spatial gra- 
dients in the temperature and the chemical potential. If 
one chooses the frame in the direction of the overall lo- 
cal energy flux, then = is concluded without the 
loss of generality. This frame is called the energy frame 
or Landau frame. From now on I consider this energy 
frame for the formulation of the hydrodynamic scheme. 

In (l+l)-dimcnsional dissipative hydrodynamics, the 
flow is expressed with the flow rapidity Yf as w M = 
(coshY/, 0, 0, sinhF/), where z-axis is the direction of 
the expansion. It follows from the orthogonality rela- 
tion tt^u^ — and the traceless condition ir^ = that 
the shear stress tensor can be expressed with a single in- 
dependent variable ir = ir 00 — tt 33 , which can be called 
shear pressure, as 

/ -sinh 2 17 -coshF/smli!) \ 
ilv I 

\ 

V - cosh Yf sinh Y f - cosh 2 Y f ) 

(3) 

Likewise, because of the orthogonality relation V^u^ = 
0, the baryon dissipation current can be expressed as 

W 1 = (-sinhY/, 0,0, -cosh Y f )V, (4) 

where V can be called baryon dissipation. 
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The energy momentum and the net baryon number 
conservations then read 



De = 

(e + Po + n - n)DY f = 



-(e + Po+n-7r)VYf,(5) 
-V(P + n-7r), (6) 



and 



DriBo = -nBoVYf + W + VDY f , 



(7) 



where the time- and the space-like derivatives in this ge- 
ometry are 

D = cosh(Yf — r) s )d T + — sinh(Y} — f] s )d„ B , (8) 

T 

V = sinh(Y) - rj s )d T + - cosh(Ff - ^ct ■ (9) 

T 

Here (r, ry s ) is the relativistic coordinate defined as t = 
t cosh r/s and z = t sinh rj s . t is the proper time and rj a 
the space-time rapidity. 

One further needs three constitutive equations to de- 
termine the space-time evolution of the systems besides 
the equation of state Po = Po(eoj^so)- Here I introduce 
the full second order dissipative hydrodynamic equations 
from Ref. [13 which extends the Israel- Stewart theory 
28] for the systems with particle number changing pro- 
cesses. The constitutive equations for the bulk pressure 
n, the baryon dissipation V, and the shear pressure tt in 
the (l+l)-dimensional geometry can be expressed as 



satisfy Onsager reciprocal relations [38[, are explicitly 
present in the equations. The diagonal baryon charge 
conductivity is related with the baryon diffusion coef- 
ficient as Db — {d[iB / dnBo)\TKv /T . The cross con- 
ductivity causes the Soret effect, which is the chemical 
diffusion induced by the thermal gradient and flow ac- 
celeration. It is related to the thermo-diffusion coeffi- 
cient as D T = (dfiB /dT) Po Kv /n B o + (kvw — ^bhv)/T. 
The ratio fcy = Dt/Db is called Soret coefficient. The 
cross coefficient can either be positive or negative while 
the semi-positive definite condition of the transport co- 
efficient matrix n V w < nv K w is satisfied. Here K\y is 
the thermal conductivity but does not explicitly appear 
in the formalism because of the frame choice, tti, tv 
and Ty are the relaxation times, xnn, Xmi an d Xvv 
the second order self-coupling coefficients, and xn-m XirV, 
and xvn the second order cross coefficients. The terms 
are combined using the conservation laws and the Gibbs- 
Duhem relation with the truncation to the second order. 
The couplings among different dissipative quantities at 
the second order can in principle have non-trivial effects 
when there are quantitative hierarchies among the dissi- 
pative quantities. 

The equations involve many time-like derivatives, po- 
tentially increasing the numerical difficulties. Here they 
are solved with an advanced version of the multiple it- 
eration algorithm I have developed for Ref. [34[ coupled 
with the piece- wise parabolic method [39| . 



£>n = i- ( - n - (nn^VY f - ( U seD^ + Cro„ B 



T 



T 



XnnnVI7 - XnvW^- - Xnv^V^ 



xZvVDY f - xEyW + xn^VYf 



DV 



1 

TV 



V + K VB W f -^- - k VjW (s7 1 



T 



A B - 1 



X^DY f + x^Vtt + XvnllV 



Mb 
T 



x£ n nv- - xvnKDYf + x^vn 



D 



(10) 



T T 



-DY t 



XvvVVYj + Xv^^r + Xy^V- 



(11) 



III. THE MODEL 

The relativistic hydrodynamic model describes macro- 
scopic motion of the fluid with the conservation laws and 
the constitutive equations. It means that the QCD equa- 
tion of state and the QCD transport coefficients must be 
given as parameter input to perform hydrodynamic cal- 
culations for the quark matter. Also the hydrodynamic 
description works in the intermediate stage of the colli- 
sions around r ~ 1-10 fm/c. Thus one further needs to 
introduce initial conditions and freezeout to link hydro- 
dynamic analyses with experimental observables. Here 
the color glass theory is employed to estimate the ini- 
tial conditions for the energy density and the net baryon 
density profiles. The hydrodynamic flow is converted to 
particles spectra via freezeout. 



Dn = — - 7T + -T)VYf 

+ X.^VYj + x.nnVF/ - x.vVV^- 

- X*vV^\ + X°vVDY f - xZvVv) , (12) 

where rj is the shear viscosity, Cnn, Cn<5e, and CnSn B 
the bulk viscosities, Ky the baryon charge conductivity 
and Kvw the baryon-heat cross conductivity. It should 
be noted that the linear cross terms, which implicitly 



A. Equation of state and transport coefficients 

The equation of state (EoS) and the transport coeffi- 
cients are the static and the dynamical responses of a 
thermodynamic system. They depend on microscopic 
properties of the medium, and are necessary input for 
hydrodynamic models. Obtaining these quantities from 
the first principle calculations is, however, generally a 
very non-trivial issue, especially for finite baryon chemi- 
cal potential systems due to the fermion-sign problem of 
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the lattice QCD calculations. Here the finite density EoS 
is constructed with the Taylor expansion method up to 
the second order, 



P(7> B ) _ P(T) 



2^4 



rp4 



2 \T 



(13) 



where \b IS the quadratic fluctuation of the baryon 
number. The latest (2+l)-flavor lattice QCD results of 
the continuum extrapolations for the EoS at vanishing 
chemical potentials [40| and for the quadratic baryon 
fluctuation [4l| are employed. It should be noted that 
the fugacity exponent hb/T has to be small for the ex- 
pansion to be formally valid. This is well motivated for 
the high-energy heavy ion collisions where the net baryon 
density is relatively small. 

The transport coefficients of the hot matter are more 
difficult to obtain from the first principle calculations. 
To the best knowledge of the author, no conclusive re- 
sults for these coefficients are available so far. This 
makes constraining the coefficients from experimental 
data in a dissipative hydrodynamic modeling one of the 
goals of heavy ion physics. Here I choose model coef- 
ficients for demonstrative purposes. The shear viscous 
coefficient is introduced from the conjectured minimum 
boundary rj/s = 1/Air [13] from Anti-de Sitter/conformal 
field theory (AdS/CFT) correspondences where s de- 
notes the entropy density. The three bulk viscous co- 
efficients Cnn, CiWe, and CuSn B are naturally expected in 
the linear response theory. £rm is the diagonal compo- 
nent in the transport coefficient matrix and the others 
are off-diagonal, or cross, components. The ratios C/ 7 ! 
are obtained from the non-equilibrium statistical opera- 
tor method with a 4 -theory estimation 43]. Following 
the discussion in Appendix [XJ I combine them to effec- 
tive bulk viscous coefficient £ = 5(| — c 2 s )rj at the first or- 
der. The coefficient phenomenologically exhibit the peak 
structure around the crossover temperature T c . 

The baryon charge conductivity and the baryon-heat 
cross coefficient characterize the finite density non- 
equilibrium processes. Using the charge conductivity 
D B = 1/2ttT in an AdS/CFT framework 0, the for- 
mer is estimated as 



KV 



Cy / d[li 



2tt \driB 



c-v 



2tt : 



(14) 



where the Taylor expansion of the EoS (1131) is utilized, cy 
is a dimensionless constant introduced for parametriza- 
tion. Here cy — 1 is considered unless specified oth- 
erwise. The cross coefficient is parametrized from the 
dimensional analyses, the matter-antimatter symmetry 
V^{hb) = — ^b), and the implication from the 
semi-positive definiteness as nyw = cvw[n-BoT/(eo + 
Po)] y/nyKw- The semi-positive definite condition of the 
transport coefficient matrix in this case is explicitly ex- 
pressed as Cy W [nBoT/(e + Pq)] 2 < 1- Here the thermal 
conductivity is chosen as Kw = 5r)T (43|. Note that the 
charge conductivity in this model remains finite in the 
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FIG. 1. (Color online) The dimensionless ratios of the shear 
viscosity (solid line), the effective bulk viscosity (dotted line), 
the baryon charge conductivity multiplied by temperature 
(dashed line), and the baryon-heat cross conductivity (dash- 
dotted line) to the entropy density at \ib = 0.05 GeV. The 
temperature region relevant to the numerical simulation 0.1 
GeV < T < 0.5 GeV is shown. 



limit of vanishing chemical potential, but does not in- 
duce the charge dissipation current out of global chemi- 
cal equilibrium because the cross coefficients vanish in the 
limit and it forms an isolated partial matrix in the full 
transport coefficient matrix. The lack of cross coefficients 
at the vanishing limit of the corresponding chemical po- 
tentials is important because it ensures that hidden con- 
serving quantities - strangeness or yet unknown charge - 
do not affect the physics. The dimensionless ratios of the 
linear transport coefficients r)/s, C,/ s, KyT/s, and Kyw/ 8 
with cy = 1 and cwv = 5 for a constant chemical po- 
tential [lb — 0.05 GeV are plotted in Fig.[TJ One can see 
that the baryon-related transport coefficients are smaller 
than the others in the current parameter settings. 

The relaxation times are also estimated from the string 
theoretical framework as tv = (2 — ln2)/27rT, rn = 18 — 
(9 In 3- v / 37r)/24 7 rT, and ry = ln2/27rT [44]. The other 
second order transport coefficients are parametrized as 
Xab = c AB [n B0 T/(e + P )] n T m T A where A and B Ac- 
note the types of dissipative currents and cab is a dimen- 
sionless constant, n = 1 is employed for the baryon- non- 
baryon cross coefficients and n — for the others, fol- 
lowing the matter-antimatter symmetry arguments men- 
tioned earlier, m is the dimension parameter chosen so 
that the dimensions of the terms in each constitutive 
equation are matched. For the most part in this study, 
however, cab = is employed to observe the qualita- 
tive nature of the dissipative processes. Note that, in the 
present analyses, the focus is on the net baryon density 
and the interplay of different dissipative currents, and it 
is beyond the scope of the paper to precisely obtain the 
transport coefficients. 
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B. Initial conditions 

The initial condition for the energy density is con- 
structed in the color glass theory. The theory describes 
the colliding two nuclei as saturated gluons called color 
glass condensate (CGC). Since the gluons are dominant 
against the valence quarks, the initial energy density of 
the hot matter is estimated from the gluon distribu- 
tion by assuming the profile of the energy distribution 
is not significantly modified during the early thermaliza- 
tion stage. The Nara Monte-Carlo adaptation [45l448l | 
of Kharzeev-Levin-Nardi model (MC-KLN) [49| is em- 
ployed for the estimation. The rapidity distribution of 
the transverse gluon energy density can be expressed in 
/cT-factorization as 



dE T 



4nN r a. 



d 2 p T 



cPrTdy iV 2 — 1 J pt 

x ip 2 (x2, (pr - k T ) 2 ), 



d 2 kT<pi(xi, k T ) 



(15) 



where N c — 3 is the number of colors, as the QCD 
coupling, X\ t 2 — Pt exp(±y)/Y / s, and cpx t 2 the unintc- 
grated gluon distributions. The lower index T denotes 
the transverse component of given position or momen- 
tum. The gluon distribution is saturated at the scale 
QIa( x ^±) = Qloi T A(r±)/T Afi }(x /x) x where A = 0.28 
is experimentally motivated. Ta is the thickness function 
to account for the nucleus geometry. Here the parame- 
ters are chosen as xq = 0.01, Ta.q = 1.53 fm -2 , and 
Ql — 2 GeV [46]. Since the hydrodynamic medium is 
expected to be constituted by low p? partons, the pt 
window 0.1-3.0 GeV is set. The space-time rapidity rj s 
is matched with the momentum rapidity y to obtain the 
initial condition in the configuration space. 

The initial condition for the net baryon density is also 
constructed from the color glass picture by assuming it 
is proportional to the valence quark parton distribution 
function [181 ]. The net baryon distribution for a nucleus 
reads 



dN 



B-B 



C 



dy 



(27T) 



d 2 p T 
Pt 



xiq v (xi)(p(x 2 ,pr), (16) 



where q v is the valence quark distribution. Here the 
normalization C is determined so that the integrated 
N B _g matches the number of participants. The dis- 
tributions are determined with the same settings as the 
no-fragmentation case in Ref. |l8[ but with A = 0.28 and 
the next-to-next-to-leading order fit results of the valence 
quark distribution (50j . Note that the net baryon distri- 
bution is sensitive to the parameters, and here they are 
chosen to yield much steeper distribution at the initial 
stage so the distribution after the hydrodynamic evolu- 
tion roughly reproduces experimental data. The tail con- 
tribution beyond the beam rapidity is exponentially cut 
off since hydrodynamic description do not apply to the 
region beyond the freezeout and the Taylor expansion- 
based EoS is not expected to work for the dense and cold 
matter. This leads to the reduction of the total baryon 



number of participant nucleons by 8.7% at the RHIC and 
0.1% at the LHC for the most central 0-5% events. 

The initial conditions for the dissipative currents are 
not well known and this is also a quite interesting issue 
by itself. Here they are chosen as non-existent at the 
initial time, i.e., n(-r ,?7 s ) = V(T ,r] s ) = 7t(t ,?7 s ) = 
for the clear view of the non-equilibrium effects on the 
fluids and also for avoiding possible overestimation of the 
viscous and the dissipative effects from ambiguity. Here 
tq is the initial time set to tq = 1 fm/c. 



C. Freezeout 

As the system cools down with the time evolution and 
becomes dilute enough, the hydrodynamic simulation has 
to be stopped and the flow field needs to be converted 
into particles. One conventionally employs the Cooper- 
Frye formula [5l[ at a freezeout hypersurface S. The 
formula reads 



9i 



d 2 N, 



d?pxdy (27r) 



(17) 



for the particle species i, where Qi is the degeneracy, da^ 
the freezeout hypersurface element, and fi the phase- 
space distribution. £ is taken as an isothermal surface 
because the chemical potential dependence of the bound- 
ary is sufficiently small for the high-energy collisions. 
This gives rise to the concept of the freezeout temper- 
ature Tf. fi can be separated into the equilibrium distri- 
bution f9 and the distortion of distribution Sfi, where 



exp 



h^B 



T 



(18) 



Here hi is the baryon number, i.e., bi = +1 for baryons, 
— 1 for anti-baryons, and for mesons, e, denotes the 
quantum statistics as ej = +1 for fermions and —1 for 
bosons. The effects of the off-equilibrium distribution 
Sfi is not treated here because its correction on the 
PT-integrated net baryon rapidity distribution would be 
small while the off-equilibrium expansion is applicable, 
since the stability condition for the net baryon density 
requires Sub = and | <^Vg | / 1 iVg 1 = V/ubo would be 
generally very small. 

It should be emphasized here that the focus of this 
paper is to estimate the finite-density hydrodynamic ef- 
fects with non-equilibrium processes on the hot medium 
because it would have the dominant effects on the net 
baryon distribution. The hadronic cascade at the later 
stage is not considered because the modification on the 
net baryon distribution during the hadronic stage would 
be small as the baryon number does not change in the 
hadronic decay and diffusion process in a hadronic gas is 
expected to be slow pjl, [HJ ■ 

The information on the chemical freezeout is implic- 
itly contained in the lattice QCD EoS in the hadronic 
phase since the hydrodynamic flow does not specify its 
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contents. If the EoS for each hadronic component was 
known, one could in principle incorporate the chemical 
freezeout explicitly by introducing conservation laws for 
the hadrons instead of the one for the baryon charge [53[ . 



IV. RESULTS 

The most central 0-5% events are considered for the 
initial conditions. The mean numbers of participants are 
357 for the RHIC and 385 for the LHC. The initial tem- 
peratures and the chemical potentials at mid-rapidity are 
T = 419 MeV and fi B = 20.6 MeV for the RHIC and 
T = 490 MeV and [i B = 6.5 MeV for the LHC. I employ 
two freezeout temperatures Tf = 0.16 GeV and Tf = 0.14 
GeV. The former early freezeout scenario is motivated 
by the possible break-down of hydrodynamic applicabil- 
ity due to the bulk viscous effects [54|. The net baryon 
distributions are constructed by taking into account the 
contributions of the hadron resonance [55| up to 2.5 GeV 
at freezeout. 



A. Net baryon distributions at the RHIC 

The net baryon distributions of the Au-Au collisions 
at ^/saT/v = 200 GeV with and without non-equilibrium 
corrections are shown in Fig. [5] for the freezeout tem- 
peratures Tf = 0.16 GeV and T f = 0.14 GeV. The off- 
equilibrium parameters are set as cy = 1, cyw = 0, 
and cab = for the moment. One can see that the hy- 
drodynamic flow tends to carry the net baryon density 
to forward rapidity, broadening the flat region with rel- 
atively small baryon density around mid-rapidity. The 
viscous hydrodynamic results are less flattened because 
the shear and the bulk viscosities prevent the expansion 
by effectively reducing the longitudinal pressure. The dis- 
sipative hydrodynamic results, which include the baryon 
dissipation in addition to the viscosities, differ from the 
viscous hydrodynamic results. The gradient in the chem- 
ical potential induces the baryon dissipation current into 
mid-rapidity region, further steepening the mid-rapidity 
valley. The data points are scaled from the net pro- 
ton distribution of the BRAHM experiments Q, which 
show relatively good agreement with the hydrodynamic 
results. It should be noted that quasi-quantitative dis- 
cussion is allowed here because the initial distribution is 
sensitive to the parameters in the color glass theory and 
also the transverse dynamics is integrated out. 

The magnitude of collisional transparency is quantified 
by the baryon stopping. The rapidity loss is defined as 
(Sy) = y p — (y) where y p is the rapidity of incoming 
projectile and 



(a) T= 0.1 6 GeV 



(y) 



y " dN 

y- 



B-B 



(y) 



dy 



dy , 



y* dN B _ s (y) 
dy 



dy. (19) 
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evolutions lead to the reduction in the average rapidity 



FIG. 2. (Color online) The initial net baryon distributions 
based on the color glass theory (solid curve) and the ones 
with the ideal (dotted curve), the viscous (dash-dotted curve), 
and the dissipative (dashed curve) hydrodynamic effects at 
(a) T f = 0.16 GeV and (b) Tf = 0.14 GeV for the Au-Au 
collisions at ^snn = 200 GeV. The experimental data points 
are the scaled results of the net proton distribution from the 
BRAHMS collaboration 0. 



loss to (Sy) = 2.09, 2.16, and 2.26 for T f = 0.16 GeV and 
(Sy) = 1.99, 2.06, and 2.19 for T f = 0.14 GeV, respec- 
tively, when that of the initial net baryon distribution is 
(Sy) = 2.67. The fact that hydrodynamic evolutions visi- 
bly reduce the rapidity loss suggests that the baryon stop- 
ping in the RHIC would deviate less significantly from 
the linear extrapolation of the AGS and the SPS results 
at the formation of the hot medium, but the net baryon 
is carried to the forward rapidity by the hydrodynamic 
medium interaction, effectively enhancing the observed 
transparency. It also indicates that the kinetic energy 
loss for the QGP production is larger, and part of the 
energy is transferred back to the net baryon component 
from the produced medium afterward. 
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The fact that the effects of the baryon dissipation cur- 
rent could be visible on the net baryon distribution is of 
importance because it suggests that one would have to 
take the diffusion process into account to quantitatively 
understand the experimental data. It would also play an 
important role in constraining the yet unknown initial 
condition of the net baryon distribution. The actual ef- 
fect of the baryon diffusion could be larger because the 
current baryon charge conductivity is moderate as shown 
in Fig. 1. It is note- worthy that though the shear and 
bulk viscosity and the baryon dissipation seem to have 
similar effects on the net baryon distribution, the former 
enhances the baryon and the anti-baryon distributions in- 
dividually while the latter only increases their difference, 
suggesting that its contribution to the averaged distribu- 
tions is small. 

Comparing the two freezeout temperatures, the broad- 
ening effect is larger for the late freezeout case because of 
the longer hydrodynamic evolution. This indicates that 
the late freezeout tends to allow larger dissipative and 
viscous coefficients. It should be noted that the total 
baryon number is slightly smaller than the initial nuclei 
at Tf — 0.14 TeV because of the Cooper-Frye formulation 
of the freezeout. The equation of state and the baryon 
fluctuation of the kinetic theory and those of the lattice 
QCD have to be identical at freezeout to perfectly con- 
serve the energy and the net baryon because they are to 
be reproduced in relativistic kinetic theory from the flow, 
the temperature, and the chemical potential. This would 
make it slightly more difficult to distinguish the freeze- 
out temperature dependence from the slopes of the net 
baryon distribution, but of course leaves the freezeout 
temperature dependence of the average rapidity losses 
unaffected. 



B. Effects of cross terms 

I next explore the effects of the cross terms in the 
equations of motion for the dissipative currents (|10[) - 
(Tl"2"j) . Fig. [3] shows the Soret effect in the QGP medium, 
which is induced by the linear thermo-diffusion term in 
the baryon dissipation, at the RHIC for Tf = 0.16 GeV. 
cyw — 5,0, and —5 are employed for the cross coef- 
ficient. The semi-positive definite condition is checked 
to be satisfied throughout the time evolution. One can 
see that the positive and the negative cross coefficients 
qualitatively lead to reduction and enhancement of the 
baryon diffusion effect, respectively. The results are con- 
sistent with the linear order analyses in Eq. (IA4[) that the 
positive thermo-diffusion coefficient effectively reduces 
the charge conductivity. The magnitude of the thermo- 
diffusion effect, on the other hand, is relatively small in 
the present calculations because the cross coefficient be- 
comes non-vanishing only in baryon rich region at for- 
ward rapidity in high-energy heavy ion collisions. This 
suggests that Soret effect might change the temperature 
and the chemical potential dependences of the transport 
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FIG. 3. (Color online) The net baryon distribution at the 
initial stage based on the color glass theory (solid curve) and 
the ones after dissipative hydrodynamic evolutions with the 
cross coefficients cvw = 5 (dotted curve), cvw = (dashed 
curve), and cvw = —5 (dash-dotted curve) at Tf = 0.16 GeV 
for the Au-Au collisions at y^sjvjv = 200 GeV. 



properties, but it would be effective only at forward ra- 
pidity. 

The second order cross terms could also affect the re- 
sults because of the hierarchy in the magnitude of dissi- 
pative currents mentioned in Sec. HH Numerical estima- 
tions with finite cab indicate that the bulk-shear cross 
term in the bulk pressure, and the baryon-shear and the 
baryon-bulk terms in the baryon dissipation can be rel- 
evant, assuming the magnitude of the transport coeffi- 
cients is roughly of the same order. Note that in general 
there is much ambiguity in the magnitude of the second 
order transport coefficients. The result is consistent with 
the fact that the shear pressure is larger than the bulk 
pressure, which in turn is larger than the baryon dissipa- 
tion in relativistic heavy ion collisions. 



C. Net baryon distributions at the LHC 

The prospects for the Pb-Pb collisions at ^/FJvtv = 2.76 
TeV in the LHC experiment are shown in Fig. 2) The net 
baryon distributions are still visibly modified by the hy- 
drodynamic flow, effectively increasing the nuclear trans- 
parency in the collision. On the other hand, the effects of 
viscosities and dissipation are much smaller. The finial 
average rapidity losses after the ideal, the viscous, and 
the dissipative hydrodynamic evolution are (6y) = 3.48, 
3.52, and 3.55 for T f = 0.16 GeV and (5y) = 3.44, 3.48, 
and 3.51 for Tj = 0.14 GeV, respectively. The initial 
rapidity loss in the current parameter settings is (Sy) — 
3.88. Comparing the early and the late freezeout cases, 
the latter tends to leave more room for viscosities and 
dissipation as is found in the RHIC settings. 
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due to the lack of explicit chemical potential dependence 
of the transport coefficients including the shear and the 
bulk viscosities in the current modeling. 

Comparing the results with the ones at the RHIC, 
one can also find that the overall hydrodynamic effect 
is smaller at the LHC. One of the reasons would be the 
fact that the mean rapidity loss in the initial net baryon 
distribution for the current parameter settings is large, 
i.e., the peak of the distribution is around y ~ 4-5 while 
the beam rapidity is y ~ 8. The pressure is larger at 
the LHC, but the pressure gradient, which drives the 
net baryon current, is relatively small in this rapidity re- 
gion. Thus the actual hydrodynamic effect could be large 
enough to be measured in the LHC experiments for the 
initial distributions with smaller mean rapidity loss or 
larger pressure gradient. 



(b)T= 0.14 GeV 



initial 

ideal hydro 

viscous hydro 

dissipative hydro 




7 8 9 



FIG. 4. (Color online) The initial net baryon distributions 
based on the color glass theory (solid curve) and the ones 
with the ideal (dotted curve), the viscous (dash-dotted curve), 
and the dissipative (dashed curve) hydrodynamic effects at 
(a) Tf = 0.16 GeV and (b) T f = 0.14 GeV for the Pb-Pb 
collisions at ^snn = 2.76 GeV. 



The smaller dissipative effect could be understood as 
a result of the smaller spatial gradients in the fugacity 
exponent /i^/T, the thermodynamic force to the baryon 
dissipation, at the LHC. The shear and the bulk viscous 
effects are reduced for a different reason, because they 
do not directly respond to the difference in the fugac- 
ity exponent. The hydrostatic pressure Po increases with 
the collision energy more than the shear and the bulk 
pressures do, and the effects of viscous corrections are 
decreased in the effective pressure P = Pq + H — 7T, re- 
ducing the difference between the ideal and the viscous 
hydrodynamic results. It should be noted that the effect 
of baryon dissipation could be larger since the charge con- 
ductivity employed in the estimations is very moderate 
as mentioned earlier. Also the off-equilibrium corrections 
on the net baryon distribution might be underestimated 



V. DISCUSSION AND CONCLUSIONS 

A relativistic dissipative hydrodynamic model of the 
high-energy heavy ion collisions at finite density which 
takes account of shear viscosity, bulk viscosity, and 
baryon dissipation with evolving flow is developed for 
the first time. The hydrodynamic framework is employed 
from the second order theory extended for the systems 
with particle number changing processes. The initial con- 
ditions for the energy and the net baryon distributions 
are constructed from the color glass theory, and the EoS 
is employed from the Taylor expansion approach of the 
finite density lattice QCD to improve quantitative accu- 
racy. I find that the average rapidity loss for the baryon 
stopping is reduced during the hydrodynamic evolution, 
which would mean that the observed transparency of 
the collision at the RHIC is effectively enhanced in the 
medium interaction. This suggests that more energy is 
available for the production of the hot medium at the ini- 
tial stage than was implied from the experimental data, 
and the strongly-coupled medium re-distribute part of 
the energy back to the net baryon components. 

The net baryon distribution could also be sensitive to 
baryon dissipation as much as to viscosities. It should be 
noted that the current dissipative coefficients are rather 
moderate as shown in Fig. [T] and the actual diffusion 
effects can be larger. The effect of baryon dissipation 
could be important for explaining the experimental data 
and constraining the initial conditions, which are not 
well known. One would need to introduce other observ- 
ables such as the transverse momentum spectra of the 
net baryon to constrain the transport coefficients of the 
hot QCD matter at finite baryon density from the col- 
lider experiments because of the ambiguities in the choice 
of initial conditions. The effects of the cross coefficients 
are also numerically investigated, and it is found that the 
thermo-diffusion effect, or the Soret effect, might mod- 
ify the magnitude of the baryon diffusion, but the ef- 
fect is limited to the forward rapidity region because the 
baryon-heat cross conductivity vanishes for the baryon- 
free medium due to the matter- antimatter symmetry. 
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The second-order cross terms among the dissipative cur- 
rents of different magnitudes would also be important for 
the quantitative analyses. The late freezeout is found to 
allow larger viscosities and baryon dissipation since the 
longer hydrodynamic evolution widens the mid-rapidity 
valley. 

A possible source of overestimation for the hydrody- 
namic effects would be the lack of transverse expansion. 
The temperature and the chemical potential tend to be 
larger for the longitudinal geometry as the energy and 
the net baryon densities cannot spread into the trans- 
verse directions. It should be noted that while the trans- 
verse dependence of the net baryon distributions is exper- 
imentally implied to be small [35( , the accelerated cooling 
would lead to the reduction in the effect of hydrodynamic 
evolution. This is partially taken into account by employ- 
ing the early freezeout scenario. The off-equilibrium cor- 
rections with transverse expansion would be more non- 
trivial and worth investigating. Also the parametrization 
of the transport coefficients for finite density systems, 
especially the conductivities Ky and KyWi needs to be 
improved through theoretical and experimental develop- 
ments for more quantitative discussion. 
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Appendix A: TRANSPORT COEFFICIENTS AT 
THE LINEAR ORDER 

The determination of the transport coefficients for a 
given system has been one of the long standing issues 
even at the linear order. The bulk pressure at the first 
order is expressed as 



n 



-C: 



f . 



mi 



T 
Cnn 
T 



T 



T 

0(S 2 ) 



dP 
de 



nir, 



T 

( dP 



T \dn B c 



= -CV M u" + 0(8 2 ). 



(Al) 



A naive (/) 4 -theory analysis [43j yields (use = — 3Cnn 
along with Cnn/T = 5i]/3. If the system is free of con- 
served charge currents, the energy-momentum conserva- 



tion and the Gibbs-Duhem relation yield 



1 



C = 5( -~c 2 ]?/• 



(A2) 



which satisfies the conjectured lower boundary £ > 2(-| — 
c 2 )ri in the M = 2* gauge theory Since the squared 
sound velocity is around I /3 except for the crossover re- 
gions, the cancellation of the linear terms lead to small 
bulk viscosity. It implies that the existence of the On- 
sager cross term is a reason for the general smallness 
of bulk viscous coefficient in hydrodynamic systems. In 
other words, the diagonal bulk viscosity Cnn prevents the 
expansion of a system but the energy cross coefficient 
(nae encourages the expansion in the effort of decreas- 
ing the temperature, leading to the overall cancellation 
of the effects. 

When the system has finite chemical potential, the 
expression is subject to non-trivial contribution from 
the density cross coefficient. For a special case where 
CnSn B = [«Bo/(eo + Po)] x Cn«5e, however, the expression 
is dramatically simplified and one again obtains Eq. (|A2I) , 
because the sound velocity is expressed as 



dP 

de 

dPo\ 
deo], 



( dP 



eo + Pq V driBo 



(A3) 



The sign of CnSn B Aip s f° r antimatter systems, correctly 
capturing the fact that the overall bulk viscosity C is sym- 
metric under the charge conjugation. I employ this form 
of bulk viscous coefficient since the main focus is on the 
net baryon density and its dissipation. 

The baryon dissipation current can be expressed in a 
similar fashion using the conservation laws and the ther- 
modynamic relation as 



K.VW 



T T 



= KV 



eo 



-K.VW 



W^ + 0{6 2 ) 



^nV^+0{8 2 ), 



(A4) 



where k is the effective baryon charge conductivity. An 
extreme case is Ky — [«so/( e o + Po)] x K vw> where the 
current vanishes altogether at this order because of the 
Soret effect. It should be noted that the vector cross 
terms arising from multiple conserved currents cannot 
be integrated in this method. In the paper, the two coef- 
ficients are treated separately to better illustrate the role 
of each linear term. 
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